Dynamics of localized spins coupled to 
the conduction electrons with charge/spin currents 
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The effects of the charge/spin currents of conduction electrons on the dynamics of the localized 
spins are studied in terms of the perturbation in the exchange coupling Jk between them. The 
equations of motion for the localized spins are derived exactly up to 0(J%), and the equations for 
the two-spin system is solved numerically. It is found that the dynamics depends sensitively upon the 
relative magnitude of the charge and spin currents, i.e., it shows steady state, periodic motion, and 
even chaotic behavior. Extension to the multi-spin system and its implications including possible 
"spin current detector" are also discussed. 
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PACS numbers: 72.25.-b, 72.25.Hg, 73.23.-b, 85.75.-d 

The interplay between the current and spins is the sub- 
ject of recent intensive studies HEHI1HI1HSH 
EE IH El HE HI- The giant magneto-resistive effect 
(GMR) 1], the spin valve effect la, the current-driven 
domain wall motion [E & IE IE 011 IE EE El El El, 
and the colossal magneto-resistance (CMR) are the 
representative phenomena due to the interaction between 
the spin state and the charge dynamics. Especially the 
charge current in the ferromagnetic system is necessar- 
ily spin polarized, and the spin torque associated with it 
drives the nontrivial phenomena [E El ■ O n the other 
hand, the pure spin current without the charge current is 
highli ghte d in the studies of spintronics (such as spin Hall 
effect Il5l 1 1(1 and quantum spin pumping (such as spin 
battery [17|). Spin current is present even in the insulat- 
ing magnets. In the latter case, the electric polarization 
is induced by the spin current associated with the non- 
collinear spin configuration, offering a new mechanism of 
the magneto-electric effect [fjj ]. 

However, the direct observation of the spin current is 
difficult because it does not appear in the Maxwell equa- 
tions. One possible way is to measure the electric field or 
voltage drop accompanied by the spin current due to the 
Aharanov-Casher effect |fjj,|20|. However, the magnitude 
of the voltage drop is estimated to be very small. There- 
fore a more sensitive method to detect the spin current 
is desirable, which we propose below. 

In this paper, we study the dynamics of the local 
spins weakly coupled to the conduction electrons with 
charge/spin currents. The equation of motion for the 
localized spins are derived up to the second order in 
the Kondo coupling Jk between the conduction electron 
spins and localized spins. This equation of motion is 
solved numerically and variety of dynamics is found in- 



cluding the steady state, periodic motion, and even the 
chaotic motion in the phase diagram of charge current j c 
and spin current . It is not unrealistic to control both 
of them by several methods proposed 0, in addition 
to the conventional spin-polarized current injection from 
the ferromagnets. It is also noted that the polarization 
of the spin current component is controllable only in the 
case of weak Jk- In the strong coupling case, the spin 
of the conduction electron is forced to be parallel to the 
localized spin at each site, and we cannot arbitrarily ma- 
nipulate its direction. 

Model and equation of motion- We start with the double- 
exchange model of continuum version treated in the 
closed time-path Green's function (CTPGF) formalism 
|2l| . which is recently applied to the spin dynamics in a 
Josephson junction. (We shall take the units in which 
h = c=l.) 



Z[S(r,t)} = Tv[T p exp(iS)}, 
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- Ut] + v JEllLs{r,t)-*{r,t) 
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dt+ - 



dt. 



(2) 
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where Z[S(r,t)] is the partition function in the CTPGF 
formalism, T p is the time-ordering with respect to the 
closed time-path, Swzwn is the Wess-Zumino-Witten- 
Novikov (WZWN) term describing the Berry phase of 
the localized spins, Ho is the Hamiltonian for free elec- 
trons. For simplicity, we consider Hq of the quadratic 
dispersion and isotropic effective mass m*. <S ex t repre- 
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sents the coupling between the system with the external 
electromagnetic field and the driving force for the spin 
current of the conduction electrons. The effective ac- 
tion for S(r,t) is given by S c s = —i\nZ[S(r,t)]. Note 
that the constraint |S(r,t)| = S = const, is implicitly 
imposed. 

Here we consider the weak coupling regime, \Jk\ <C 
cf (cf : the Fermi energy) and retain the terms up to 
the quadratic order in Jk- Then, as in Ref.[22T]. S c g is 
estimated as follows. 

Scs = Swzwn+SMb / dt drB(r,t) ■ S A (r,t) 

J — oc J 



where / (e&) is the distribution function in equilibrium, r 
is the electron lifetime, Vf. is the group velocity, E is the 
external electric field, e s = 1/2 and E s>± is the (virtual) 
driving force for the spin current, i.e., 



[3 >Jk=0 
— J.3 



-E, (f 



>Jk=0 



-E s 



where 



n e = k F /(3ir) 2 and kp is the Fermi momentum. 

Under the above assumption the equations of motion of 
discrete spins are obtained by using the relation S(r) <-> 
T,n S ( r -r n )S rn . 



dSi 
dt 



dt 



dtdt' / drdr' 



VucJk 



xS£(r,t)Il"(r,t;r',t')S» c (r',t% 



(4) 



Si x 



r sc \' 



,,jX Sj 



where S c (r,t) = [S(r,t+) + S(r,t_)]/2, S A (r,t) = 
S(r,t+) - S(r,t_), and n*(r,t;r',f) = -*0(t - 
t'){[a^{r,t) 1 (j v {r' 1 t l ))j K= o is the retarded spin Green's 
function of the conduction electrons in nonequilibrium, 
which represents the response of the spin polarization of 
the conduction electrons to the localized spins in the lin- 
ear order, i.e., 

M*->*)> = f dt' Jdr'U^(r,t;r',t')S^r',t'). 

(5) 

The variational principle 8S c s/5S A (r,t) — at 
SA(f,t) = gives the equation of motion for localized 
spins, 



S <i X S n 



dt 



-l 2 Si x 



0) 



where r, 



rj and 



S r JS. The time scale is 
normalized as t = 7 2 e^t where 7 = \/~Sv uc n e Jx/(2ep ). 
In the above equation of motion, we have retained the 
time derivatives up to n < 1. The detailed forms of di- 
mensionless functions ^occ( r ) an< ^ -^sc ( r ) f° r the n-th 



time-derivatives are given elsewhere |23j. F^°\r) is the 
Ruderman-Kittel-Kasuya-Yosida (RKKY) range func- 



tion,^ (1) (r) = -f 



1 — cos(2fc F r) rri( n ) 



dS(r,t) 
Ft 



S(r,t) x 



gfi B B{r,t) 



Vuc Jk 



<x(r,i)) 



and Ft^ir) oc 



(6) 



where we have used the abbreviation S c (r,t) — > S(r,t). 

For simplicity, we shall consider the case B(r,t) = 
hereafter. In order to estimate Hj}„(r,t;r',t'), 
we approximate the anti-commutation relation 
jjr=o by its c-number value in 
equilibrium, and the nonequilibrium nature of the 
system is taken into account through the nonequilibrium 
distribution function g(ek), where = k 2 /{2m). 

(cl(r',t')c a (r,t))j K=0 



dk 

W) 1 



e^-^-^-^-^g^iek). (7) 



Here we adopt the following form of g{tk) in order to 
express the state with the charge current j c and the spin 
current j s of the conduction electrons. 



l(e k ) = f(e k ) + 



erv k ■ E + e s rv k ■ E s a. 



df , ; 



r ' (j s )Jk=Q- It is noted that at 
least two spins are needed in order to feci the flow of 
the charge/spin currents. This results is roughly con- 
sistent with the result obtained by Barnes for the dy- 
namical RKKY problem j^. It can be shown that the 
above equations of motion is valid in the present situ- 
ation 7 2 <C 1 and the magnetic field is not too strong 

E3 

F{f\r)-terni is the intrinsic torque between the lo- 
calized spins due to the RKKY interaction. F o (1) (0)- 
term represents the Gibert damping due to the Fermi 
bath of the conduction electrons, and F^ (r)- term 
(r 7^ 0) works as a non-local damping/puming |24l |25|. 

Fc? (r)-/Fsc (r)-term is the intrinsic torque due to the 
charge/spin currents. Depending on the separation of 
spins and the flow of charge/spin currents, F C c (r)-term 
renormalizes the original toruqe term but Fs? (r)-term 
introduces a torque term with a different symmetry, be- 
cause the spin current is not a vector but a tensor. As we 
shall see later, this term introduces the spin anisotropy 
in the direction of the spin component of the spin cur- 
rent. Fee {r)-/ Fsc\r)-teTia is the correction to the non- 
(8) local damping/pumping term due to the charge/spin cur- 
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rents. FgP (r)-term has a different symmetry from that 
of the original term and leads to a spin-anisotropic damp- 
ing/pumping. 

Dynamics of two-spin system- Now we consider the dy- 
namics of the two-spin system in detail. For the numeri- 
cal study, we normalize the charge/spin currents as 



J" 



\2 
A F 

2ee F 



[3 )j K =o, 



2e s €F 



'Jk 



,o,(10) 



where Xf — 2Tr/kp, and, take the forms of them as 
jn _ jnj Q an( j js" _ J s n %jo, where jo and n s 
are the unit vectors representing the directions of the 
charge/spin-current flow and the spin-current polariza- 
tion, respectively. The unit of charge current, 2ecp/X F , 
is estimated as ~ 10 5 Acm -2 for a semiconductor with 
n e ~ 10 18 cm -3 and cf ~ 1 meV, and as ~ 10 10 Acm -2 
for a metal with n e ~ 10 21 cm -3 and ep ~ 1 eV. The unit 
of time, 1/(j 2 cf), is estimated as 10x(l meV/e F ) ps for 
7 = 0.05. We put the two spins along j which is set to 
be the x-direction. n, is set to be the z-direction. 




FIG. 1: Phase diagram of the dynamics of the two-spin sys- 
tem in the plane of charge current J™ and spin current J s . 

Fig. H show the " phase diagram" for the dynamics 
of the two-spin system with kpr = 1.0 and r y 2 = 0.05. 
Figs. 12141 show examples of motions in the regions II- IV 
respectively. The black line represent the inner product 
of the two spins Si ■ S 2 - The solid lines of green and or- 
ange are the z-components of upstream and downstream 
spins respectively, i.e., Sf and S$. The dashed lines of 
green and orange are the x-components of upstream and 
downstream spins respectively, i.e., Sf and Sf . 

In region I, the two spins reach the final static 
collinear configuration with the easy-axis anisotropy in 
the direction of n s . In the special case of J s — 0, i.e., the 
pure charge current, the static collinear configuration is 
realized with no spin anisotropy. This region contains 
the cases with any amount of the weakly spin-polarized 
charge current and a small amount of the perfelctly spin- 
polarized charge current. In region II and III, the spins 
show the periodic motions in which the angle between 
the two spins and the z-component of each spin takes a 



constant value. The weakly-charged spin current belongs 
to region II, which is our main concern. An example of 
this case is given in Fig. where the lines for Sf (green 
solid) and Sf (orange solid) are almost overlapped. 
When | J™ | is sufficiently smaller than \J S \ the angle 
6*i2 between the two spins is a increasing (decreasing) 
function of \J S \ (|J™|). The z-componet of a spin is 
a increasing (decreasing) function of \J n \ (|J S |). The 
frequency of motion increases with both of | J n \ and | J s \. 
In the case with J™ = 0, i.e., the pure spin current, the 
static coplanar configuration with the hard-axis parallel 
(or the easy-plane perpendicular) to n s is realized and 
the angle between the two spins, #12, is determined by 
the magnitude of J s as | tan(9 12 | = \F^ } (m) / F (0) (n 2 )|. 
Fig. |21 shows another type of a simple motion which 
occurs in region III. In region IV observed the chaotic 
or complicated motion in which both the angle between 
the two spins and the z-component of each spin are 
changing in motion and do not reach a steady state. 
This region contains the cases with a large amount of 
the strongly spin-polarized charge current and a large 
amount of the strongly charged spin current. Although 
we observe various types of motions in this region, any 
sharp boundaries between them are not identified at 
least in the present simulation. We shall not undertake 
their precise classification. However, a relatively clear 
change can be observed around the dotted line in Fig. ^ 
Under the dotted line in this region, the upstream spin 
slowly precesses around n s with a cycloidal trajectory 
and the downstream spin quickly precesses around the 
upstream spin. Deep in the region IV above the dotted 
line, the motion becomes chaotic as shown in Fig. 0] 
Note that a large amount of strongly spin-polarized 
charge current leads to this class of motion, where the 
strong coupling theories ( for 7 2 3> 1 ) 0, HiJ predict 
the instability of the ferromagnetic order. Some remarks 
are in order here. First the dynamics highly depends 
on the initial condition in the vicinity of the boundary 
between the regions II and IV. Secondly there is a very 
narrow region between the regions I and IV, while the 
regions I and IV appear to touch each other. In this 
narrow region, we observe the motions similar to those 
in the regions II and III. Thus, the regions II and III 
may be connected through this region. 
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FIG. 2: Example of the motion in the region II. J n = 0.1 
and J s = 2.0. For the meaning of each curves, see the text. 
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FIG. 3: Example of the motion in the region III. J n = 2.0 
and J a = 0.6. 
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FIG. 4: Example of the motion in the region IV. J™ = 1.0 
and J s = 1.0. 

Discussion- The above consideration on the case of 
pure spin current can be easily generalized to the one- 
dimensional periodic array of spins along the direction 
of jo. The equilibrium interaction i 7 ^ ' (ao) between the 
nearest neighbor spins with the separation ao is assumed 
to be ferromagnetic. Even in the case with J s ^ but 
J™ = 0, we can construct an energy functional E c g which 
should be minimized for the spin configuration of final 
steady state. 

Ees = [> (0) foi)$ ' Sj + F^im) ■ [Si x S 3 .]' . 

i<j 

(11) 

This suggests that even the dissipative transport spin 
current works on the localized spins in the same way as 
the equilibrium spin current does, while it is not the case 
for the charge current. In equilibrium, i.e., J n = J s = 0, 
the spins align ferromagnctically, and the direction of 
alignment is arbitrary because there is no spin anisotropy 
in our model. When J s is introduced, the spins feel 
Dzyaloshinsky-Moriya-like interaction in addition to the 
ferromagnetic one. For simplicity, we shall take into 
account only the nearest neighbor interaction and set 
n s to be in the positive z-direction. Then, the varia- 
tional energy including the lowest energy state is given 
by E q = -F (0) (a )cos(qa ) + i^ c 0),2 (a ) sin(ga ). The 
lowest energy state is realized by the real-space config- 
uration Si/S = [cos(go^i + 4>o)i sin(<7oJ"j + </>o),0] where 
tan(qoao) = — Fj$' z {aa)/FQ (oo) an d 0o is an arbi- 
trary constant. The spins is expected to show the he- 
lical order and the hard-axis anisotropy parallel (or the 
easy-plane anisotropy perpendicular) to n s is introduced 
by the pure spin current. The pitch angle, |go a o| ? of 



the helix increases with the magnitude of \J S \, because 
|-Fsc^(oo)| oc |J S |. Hence this helical order is realized in 
the steady state with the pure spin current, this might be 
utilized as the spin current detector. In other cases than 
discussed above, we must consider the dynamical equa- 
tions of motion directly. However, the previous study 
of the two-spin system gives us the following insights. 
When the small amount of J™ is introduced to the pure 
spin current, the spins arise from the coplanar helical 
order in the easy plane and show the periodic motion 
around n s as in the motion of the two spins in the region 
II. In the case with a small amount of the weakly spin- 
polarized current which belongs to the region I, the spins 
align ferromagnetically as in equilibrium, but in this case 
the easy-axis anisotropy parallel to n s is introduced by 
the spin-polarized current. However, for a large amount 
of the strongly spin-polarized charge current belonging 
to the region IV, the instability of ferromagnetic order 
would be observed. The details of the dynamical motion 
of multi-spin system will be reported elsewhere. 

In conclusion, we have derived the equations of mo- 
tion of localized spins weakly coupled to the conduction 
electrons with charge/spin currents. The dynamics of lo- 
calized spins are affected significantly by these charge and 
spin currents of conduction electrons, and the dynamics 
of the two-spin system has been studied in depth by nu- 
merical calculation. The phase diagram for the dynamics 
in the plan of charge and spin currents are divided into 
4 regions, where the steady state, periodic motion, and 
even chaotic motion are observed. 
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